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Recently, the interest in local lattice actions for chiral fermions has revived, with the proposition 
of new local actions in which only the minimal number of doublers appear. The trigger role of 
graphene having a minimally doubled, chirally invariant, Dirac-like excitation spectrum can not 
be neglected. The challenge is to construct an action which preserves enough symmetries to be 
useful in lattice gauge calculations. We present a new approach to obtain local lattice actions for 
fermions using a reinterpretation of the staggered lattice approach of Kogut and Susskind. This 
interpretation is based on the similarity with the staggered lattice approach in FDTD simulations 
of acoustics and electromagnetism. It allows us to construct a local action for chiral fermions which 
has all discrete symmetries and the minimal number of fermion flavors, but which is non-Hermitian 
in real space. However, we argue that this will not pose a threat to the usability of the theory. 



I. INTRODUCTION 

The notorious fermion doubling problem has always 
complicated calculations of lattice field theories in which 
strict chiral symmetry of the fermions is required. The 
no-go theorem of Nielsen and Ninomiya ensures that ev- 
ery local, unitary and chirally symmetric lattice theory 
will have a degenerate excitation spectrum, with at least 
two independent fermion flavors [U [3]. Eliminating 
these fermion doublers requires explicit breaking of chi- 
ral symmetry on the lattice, as proposed by Wilson pQ. 
Although chiral symmetry is restored in the continuum 
limit, this situation is often disadvantageous. Other solu- 
tions to circumvent the fermion doubling problem include 
overlap fermions, which are highly nonlocal, or domain 
wall fermions, which require an extra dimension. Both 
approaches are thus computationally expensive. 

A possible solution of a different kind is the use of 
the rooting procedure. If the action is quadratic in the 
fermion fields, as is the case in lattice gauge theory, they 
can be integrated out. The resulting determinant, de- 
scribing Nf fermion flavors, can be exponentiated to a 
power l/Nf. This reduces the number of fermions to one 
in all orders of perturbation theory. Although lattice 
gauge theory constitutes the majority of all lattice calcu- 
lations, this is a restricted approach, which is also criti- 
cized of introducing non-perturbative artefacts EJ [7] . 

The rooting procedure is virtually always used in com- 
bination with the staggered lattice formulation of Kogut 
and Susskind [U [9j [10] . This is a simple construction to 
reduce the number of fermion flavours in the lattice the- 
ory. However, on a four dimensional Euclidean lattice, 
not two but four fermion flavours survive the reduction. 
The observation of a minimally doubled Dirac-like ex- 
citation spectrum in graphene has triggered the search 
for new lattice actions for local, chiral fermions with the 
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bare minimum of different flavours [TT1 [HI [HI IH1 [T51 ITfj] . 
However, these all violate at least some of the discrete 
symmetries of the continuum theory, and interactions 
will generate unwanted relevant operators in the effec- 
tive continuum theory [ 17| . This paper tries to offer a 
workaround for this problem. 

Since the presence of fermion doublers is a property 
of the equations of motion, it can be investigated at the 
classical level. The fermion doubling problem is often 
attributed to the first order derivatives appearing in the 
Dirac Hamiltonian or action. However, first order partial 
differential equations also occur in acoustics and electro- 
magnetics. Although a naive discretization also results 
in an analogue of the doubling problem, this problem can 
be completely circumvented. The solution is, not inciden- 
tally, called the staggered grid approach and is commonly 
used in numerical simulations using the FDTD (finite- 
difference time-domain) method. We apply this machin- 
ery to the equations of motion of the Dirac fermion. Re- 
moving all doublers with this approach is still forbidden 
by the Nielsen-Ninomiya theorem. The reason for this 
is the spinor nature of the fields, which have no natural 
position on the lattice, in contrast with the scalar and 
vector fields in acoustics and electromagnetics. 

In addition, formulating an action on a completely dis- 
cretized spacetime lattice without introducing new dou- 
blers is a non-trivial procedure. We illustrate a possible 
approach, which results in an action with no time dou- 
blers and all required discrete symmetries. This action 
is non-Hermitian in real space, but will not introduce 
unphysical effects in the calculation of expectation val- 
ues. We start with (1 + l)-dimcnsional fermions, where 
the doubling problem only appears when both space and 
time are discretized. The Kogut-Susskind construction 
can completely eliminate the fermion doubling problem 
on the Hamiltonian lattice, so we try to discretize time 
without increasing the number of doublers. We pay spe- 
cial attention to the conservation of the discrete symme- 
tries P, C and T, and briefly comment on how to incorpo- 
rate gauge fields. The second section considers the more 
interesting case of fermions in (3 + 1) dimensions. 
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II. STAGGERED LATTICE IN (1 + 1) 
DIMENSIONS 



At first, the original Kogut-Susskind construction on 
the Hamiltoninan lattice in (1 + 1) dimensions is re- 
viewed. This is the only simple approach which escapes 
the Niclscn-Ninomiya theorem, since only one dimension 
is discretized. By reinterpreting this approach from a 
point of view borrowed from the field of FDTD simula- 
tions, we illustrate how to construct a lattice action which 
has all the required symmetries without introducing ex- 
tra time doublers. The way to circumvent the Nielsen- 
Ninomiya theorem this time, is by giving up Hermiticity. 
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FIG. 1: Unit cells for the discretization of the pressure and ve- 
locity fields in acoustics (a), and for the electric and magnetic 
field in electromagnetics, the famous Yee cell 18 (b). 



A. The Kogut-Susskind construction on the 
Hamiltonian Lattice 



The staggered fermion formulation on the Hamiltonian 
lattice in (1 + 1) dimensions starts by introducing an aux- 
iliary lattice with twice the number of sites as the physical 
lattice. On each site k, a one-component fermionic vari- 
able 4>(k) is introduced, which is related to the continuum 
fermion field by 



<f>(2k- 1) 



■4> 2 (ka x ) 



<t>(2k) 



(1) 



The constant a x will always denote the lattice spacing of 
the physical lattice. The Dirac Hamiltonian 



H 



dx %jr{x) 



-ia- h m[3 

dx 



ip{x) 



(2) 



with a — <j x and = a z , is discretized using an asym- 
metric definition for the derivatives 

d-01/, \ ipi((k+ 1)0 - ipi(ka x ) 

——{ka x ) = , (3a) 

da; a x 

^2 n \ i>2(ka x ) - i/j 2 ((k- l)a x ) 

-j — {ka x ) = . (3b) 

da; a x 

This results in a pair of equations of motion which have 
non-trivial solutions if, in Fourier space, 



B. Reinterpretation of the finite differences 

The auxiliary lattice and the asymmetric definition of 
the finite differences in ([3]), appear to be rather random. 
Furthermore, this approximations for the derivatives are 
only accurate up to order 0(a x ), while a symmetric def- 
inition would be accurate up to order 0(a x ). This defi- 
ciency can be cured by a simple reinterpretation of the 
staggered fermion construction. 

In FDTD simulations of acoustics or electromagnetics, 
which are also governed by first order partial differential 
equations, the concept of a staggered lattice designates 
that different fields are discretized at different positions. 
FIG. Q] illustrates the standard unit cell for both areas. 
Essential to this approach is the fact that the equation of 
motion of one field is not dependent on spatial derivatives 
of that same field. 

Applying the same technique to the Dirac equation, we 
immediately arrive to a lattice scheme where the ipi field 
is discretized half in between the discretization points of 
the ^2 field. Thus, the Dirac Hamiltonian ^ can be 
discretized as 



H =a x ^ip\(ka x ) 



-i^j-^(ta x ) + m^pi(ka x ) 



+v4((*+ 1 / 2 )°x) 



1^(^ + 1/2)0 

da; 



(6) 



mip 2 ((k + 1/2)0 1 



= E(k) 2 with E(k) = Am 2 + sin 2 («;/2) . (4) with derivatives which are approximated by 



There is only one low energy sector and thus only one 
fermion flavour. A chiral transformation can be imple- 
mented as 

7 5 : ipx(ka x ) ^ ip 2 (ka x ), ip2(ka x ) -> ipi((k + 1)0 > 

(5) 

which is equivalent to the continuum chiral transforma- 
tion up to order 0(a x ). More details about this construc- 
tion can be found in [S]. 



^((k + 1/2)0 = Mfc + iK^M^ , (7a) 

da; a x 

dip2,, , ^ 2 ((fc + 1/2)0 -?/; 2 ((fc- 1/2)0 , 7M 

-7—{ka x ) = . (7b) 

da; a x 

These are accurate up to order 0(a x ). The derivatives 
are now symmetric and the arbitrariness is shifted to the 
allocation of the two field components to one unit cell, 
which has no physical content and is only necessary if 
we want to label the lattice fields by integers. With the 
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assignment 



ipi(ka x ) 



Xi (fe) 



^((fc+VsK) = 



X2(fc) 



(8) 



we arrive at the same theory as the Kogut-Susskind for- 
mulation on the mathematical level. The discrete-space 
Hamiltonian of the Xi fields is given by 



k 



■ X2(k) - X2(k - 1) , ' 

-l h mxi(k) 



n., 



.Xi(fc + 1) - xi(fc) 

-i mx2\kj 

a- 



(9) 



The difference in physical interpration between (TJ and 
([8]) only comes in to play when adding extra fields such 
as gauge fields. The best solution is then to discretize 
the gauge field twice as fine as the fermion field. There 
will be two gauge units cells in each fermion unit cell. 
This is equivalent to adding a gauge variable to each link 
in the auxiliary lattice of Kogut and Susskind. Wilczek 
pointed out why a finer discretization of the gauge field 
is an appropriate concept |19j . We will comment in more 
detail on the incorporation of gauge fields for the case of 
the spacetime lattice. 

Before continuing this discretization to the dimension 
of time, we elaborate on the symmetries of the discrete- 
space Hamiltonian We focus on the discrete sym- 
metries, which are, besides chiral invariance, P (invari- 
ance under parity transformations), C (invariance under 
charge conjugation) and T (invariance under time inver- 
sions). The chiral transformation on the lattice was al- 
ready defined in ([5]). With this new interpretation, we 
see that it is not possible to simply swap the spinor com- 
ponents "01 and ip2 in the lattice model, since they are 
defined at different points in space. As a consequence, 
this swapping should be accompanied with a shift of the 
lattice over half a lattice spacing. This operation trans- 
forms tpi to ip2 in the same unit cell, but shifts ^2 to ipi 
in the next unit cell (or vice versa, depending on the di- 
rection of the shift and the definition of the unit cell) . It 



is obvious that the kinetic term in is invariant under 
this transformation, while the mass term changes sign. 
The new definition given in ([8]) thus renders the idea of a 
chiral transformation which deviates from the continuum 
prescription more acceptable. 

The parity transformation is defined in the continuum 
theory as ip( x ) ~ ¥ /3V'( — x ) an d is nicely mapped to the 
lattice variables as 

P: Xi(k) ^ Xi(-k) , X2(fc)--Xa(-*-l), (10) 

which leaves (9) invariant. Charge conjugation in (1 + 1) 
dimensions is defined as if> — > oap* , and thus also depends 
on swapping the fermion components. Using the same 
prescription as for chiral invariance, the Hamiltonian ^ 
is invariant under 

C: Xi(k)^xt(k), Xa(*0-xl(*+l). (11) 
Finally there is invariance under time reversal, which is 
an antilinear operation. In the (1 + 1) dimensional con- 
tinuum time reversal is established as ip{x) — > /3ip(x), 
while the antilinearity results in a complex conjugation 
of the prefactors. The Hamiltonian (|9| is indeed invari- 
ant under 

T: Xi(fc)->Xi(fc), xa(fc) -> -* 2 (fc) 

+ complex conjugation of prefactors. (12) 

In the next paragraph, we try to discretize time, by 
going to the action formalism, as is used in Monte Carlo 
simulations of lattice field theory. The discrete symme- 
tries P, C and T are difficult to preserve in an action 
formulation with discrete time, without introducing new 
doublers. These symmetries are however of key impor- 
tance to ensure a correct continuum limit when interac- 
tions or gauge fields are added. Absence of these sym- 
metries would allow the appearance of new relevant op- 
erators in the effective continuum theory |17j . 



C. One-fermion chiral action on the spacetime lattice 



In FDTD simulations, time is discretized as well. The time evolution of the initial configuration of the fields is 
approximated using a leap frog algorithm, which implies that the different fields are also staggered in the direction of 
time. With the correspondence 

Mka x Ja t )= ) ^±, M(k+ 1 /2)a x ,(l + 1 /2)a t ) = ^ 2 ^ 1 (13) 



the FDTD equations for the free Dirac field would be 

f h [Xi(fc, - Xi(k, l-l)]=-± [ X 2(k, l -l)- X2 (k + f [xi(k, I) + X i(k, I - 1)] 

I h [X2(k, I) - X i{k, l-l)]=-± [ X i{k +1,0- Xi(k, I)} - f {X2(k, I) + X 2(k, 1-1)}. 



(14) 



These equations of motion result in a unique fermion flavour. The objective is thus to construct a discrete spacetime 
action which returns these equations as extremization conditions. But firstly, we determine the dispersion relation of 
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FIG. 2: Dispersion relation u(k) for ma t = 1 and ^ = a t (left) and two values of ui, for « = and K = n/2, as function of ma t 
(right). The full line represents the continuum solution w = (m 2 a 2 + « 2 a 2 /a 2 ; ) :L/ ' 2 , the dashed line displays the solution (161, 
while the dotted line illustrates the effect of having a dispersion relation 1 17 1. 



these equations. After Fourier transforming x-i(k, I) = exp(ink — ujl)xi{V) K ), non-trivial solutions only exist if 

4 sin 2 (w/2) - m 2 cos 2 (w/2) - 4r sin 2 (K/2) = . (15) 



Although the w-dependent factor in the term proportional to m looks strange, the corresponding solution 



uo = ±2 arcsin 



'm 2 a 2 + 4sin 2 (ft/2)a 2 /a 2 . 
4 + m 2 af 



(16) 



is, for ma t sufficiently smaller than 1, closer to the continuum solution u>/a t — ±(m 2 + K 2 /a 2 ) 1 ' 2 than the solution of 
the expected dispersion relation 



\ sin 2 (w/2) - m 2 » sin 2 ( K /2) = . 



(17) 



as is illustrated in FIG. [2j In particular, the dispersion relation (15 1 yields real solutions of ui for all values of k in 
the domain [— tt, +tt[. 

In order to find a discrete-time action, we can try a naive discretization of the continuous-time action corresponding 
to the discrete-space Hamiltonian([9| on a staggered spacetime grid. We temporarily restrict our attention to the 
masslcss Dirac fcrmion field for brevity of the expressions, which results in the following Minkowski action 



x\(k,i) + x \(k,i-i) 



k,l 



xi(M)-xi(M-i) , . x2(M-i)-X2(fc-i,z-i) 

at a x 



xl(k,l) + xl(k,l-l) 



\ X2(k,l)- X 2(k,l-l) . X i(fc+l,f)-xi(fc,0 
a t a x 



(18) 



This action has all discrete symmetries of the continuum theory. However, due to the averaging of the xj fields, there 
will be a line of poles for to = n in the fermion propagator, which is even worse than an extra fermion flavour and 
is thus highly unphysical. Correspondingly, the extremalization conditions dS/dx\ — do not reduce to the optimal 
leap frog equations (14 1. To circumvent this problem, averaging should be avoided. An arbitrary choice for one of 
the two positions xl(k, Z — 1) or x\{k, I) results in 

S = a t J2x\(k,l-l 



k.i 



■ Xi(M)-xi(M-i) , . x 2 (M-i)-X2(fc-U-i) 



+ Xl(kj-l) 



at a x 

■ X2(fc,/)-x2(fc,t-i) + . xi(fc + i,0-xi(fc,0 

a t a x 



(19) 



5 



This time the extremalization conditions dS/dx] = 
are given by (14) as wanted. The price we have to pay 



is the loss of Hermiticity at order 0(at). However, the 
extremalization conditions dS/dx\ = and dS/dxi — 
do not conflict, since 



dS 



dS 



dXi(k,l) \dx\{k,l-l) 



(20) 



As a consequence, both extremalization conditions re- 
sult in compatible equations of motion for x and x^ > but 
shifted over one time step. The inverse fermion propaga- 
tor in momentum space is given by 



D~ 



at 

-e il 



J / 2 sin(w/2) 
' sin(«;/2) 



iK / 2 sin(«/2) 
^ 2 sin(w/2) 



(21) 



which is Hermitian if not for the global factor exp(iw/2). 
The action ( 19 1 does not only sacrifice Hermiticity, it 



requires a special treatment concerning symmetries in ad- 
dition. A chiral transformation is now accompanied by a 
shift over half the lattice spacing, both in the direction 
of space and in the direction of time 



5 . jxi(M) -> X2(M)> 

7 ' |xa(M)->xi(* + M + i)- 

A parity transformation is still given by 

p . fxi(M)->xi(-M), 
' U 2 (fc,i)->- Xa (-fc-i,z). 



(22) 



(23) 



The action (19) is invariant under both transformations. 



It is however not invariant under charge conjugation 

Xi(k, I) -> xl(k, I) , Xa(*, - x\(k + 1,1 + 1), 
or time inversion 



xi(M) ->xi(fc,-0 > X2(M) ->• -X2(fc, 

+ complex conjugation of prefactors. 



-/-I), 



Only a combined CT transformation maps the action ( 19 ) 
back to itself 



xi(M) 

X2(k,l) 



-xl(k,-l-l), 



X \(k + 1,-1-1) 



+ complex conjugation of prefactors. 

The lack of time reversal invariance can easily be seen 
from the structure on the action, which is depicted in 
FIG. [3j The time derivative of ipi is defined at the 
position half in between the points were the two fields 
■01 appearing in the derivative are defined. We have no 
tp\ at that position to complete the discrete analogue of 
V'lffV'i m the action. Averaging leads to a line of poles 
in the inverse propagator, while adding a single nearby 
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FIG. 3: Unit cell in spacetime for the discretization of the 
fermion fields on a staggered grid. 



ip\ breaks time reversal symmetry. Under the suggested 
transformation laws for charge conjugation or time rever- 
sal, the expression for S with front factor x\{k, £ — 1) is 
transformed to the equivalent expression with front fac- 
tor xl(k,l)- Another possible solution consists of adding 
extra discretization points for the fields at the intermedi- 
ate position. This results in time doublers and is exactly 
the approach we are trying to avoid. 

However, there is a way to restore the symmetries by 
allowing more deviations of order 0(a) between the lat- 
tice and the continuum theory. In all the proposed trans- 
formation laws, the transformation of the fields x\ IS lm ~ 
plicitly assumed to be the complex conjugate of the trans- 
formation of the fields Xi- By dropping this restriction, it 
is possible to define a time inversion and charge conjuga- 
tion transformation which independently map the action 
( 19 1 to itself. This is allowed, since this lattice action is 



only meaningful in the definition of a partition function, 
in which x% an d x\ are independent integration variables. 
It makes no sense to define x\ 83 the conjugate momen- 
tum of Xi an d to define a (anti)commutation relation on 
both fields. Canonical quantization is only possible in a 
Hamiltonian framework, in which the time dimension is 
continuous. Only for a t — *• will the transformation laws 
match and the Hermiticity of the action be restored. The 
correct symmetry transformations can easily be deduced 
by further implementing the idea of staggered lattices. 
If we discretize ipj in between two ipi's, the perfect dis- 
crete analogue of V'JftV'i can be constructed. In order to 
retain action (1 1 91) in this new interpretation, we define 



Mka x ,la t ) = ^, 
4(ka x ,(l + y2)a t ) = ^, 

M(k + V2K,(z + y2K) = ^ 

V4((fc+ yajox, (i + l)at) 



(24) 



Alternatively, we could have used ipl((k + 1 /2)a x ,lat) in 
the last line, but this would require a small modifica- 
tion in the second term of action (19 1, as well as a sep- 



arate transformation law for x\ i n case of a chiral trans- 
formation. Again, this renewed interpretation does not 
change the mathematical content of the model. How- 
ever, at the physical level, it explicitly destroys any hope 
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for a Hermitian action, since the complex conjugate of 
the fields can not be defined on the same lattice. In 
other words, at the physical level, the complex conju- 
gate Xi(k, I) only equals x\(k, I) up to corrections of order 
0{a). Furthermore, this point of view can result in ac- 
tual differences when adding interactions and it changes 
the way in which correlation functions should be inter- 
preted: (x\ (k, l)xi(k' , I)) is not an equal-time correlation 
function, except for a t — > 0. And last but not least, it fa- 
cilitates the search for symmetry transformations. Using 
( 24 1 , we define a time reversal transformation as 



fields 



T 



'xi(M) -> xi(k,~l) , 
x \{k,l)^ X \(k,-l-\), 
Xz(k,l) -> -X2(k,-l - 1) , 
xl{k,l)^-x\{k-l-2), 
+ complex conjugation of prefactors. 



(25) 



This is indeed an invariant transformation for the action 
in (19). To define a charge conjugation transformation, 



which maps ipi to i/j 2 and vice versa in the continuum 
limit, we have to shift the lattice over half a lattice spac- 
ing in the space direction. This implies 



Xi(k, l)->xUk,l-l), 

*t(M) -> Xa(M) , 
X2(k,l)^ X \{k + l,l), 
k x\(k, I) -» Xx (k +!,/+!), 



(26) 



which maps (19) to itself as well. 

As a conclusion, we point out that the problematic 
term in the fermion action is ip^^tp. It does not change 
shape under unitary transformations of the Dirac spinor 
i\) and is thus always present. Treating this term on a 
discretized spacetime lattice with local approximations 
for the time derivative, will require to chose between a 
number of possibilities which all have disadvantages. Av- 
eraging is the worst solution, leading to a line of poles in 
the fermion propagator. Spreading the derivative over 
two lattice sites introduces time doublers. The only way 
to circumvent this problem is by giving up on Hermitic- 
ity of the action and a strict relation between the fields 
and their complex conjugate on the lattice. This is al- 
lowed, since the complex conjugate fields are independent 
variables in a partition function. 

And finally, it is possible to show that the non- 
Hermiticity of the action does not cause any unphysical 
effects. By performing independent Fourier transforma- 
tions on the fields x% an d xIj given by 

fxi(M) = £ KW Xi(^)e i ('") , 

x\(k,l) = E K .X t i(«^)e^ i(Kfe ^ (/+1/2)) , (1?7] 
X 2 (fc,0 = E K .X2(^^)e i (^+V 2 )-(^)) j I > 

W(M) = E_x 2 (^)e- i(K(fe+1/2) -^ +1)) , 

the original action can be made Hermitian in the new 



Xi(«,w) 


t ■ 


_X2 («,<*>)_ 





-|sm(u;/2) -£sin( K /2)' 
-^sin( K /2) -£sin(w/2). 



Xi(k,uj) 

X2(K,UJ) 



(28) 



Furthermore, it is now appropriate to consider the fields 
Xi(k,u>) and xj(/c, w) as Hermitian conjugates at the 
physical level. This way, the non-Hermiticity of the ac- 
tion in real space can be regarded as a consequence of 
performing non-matching Fourier transformations on the 
fields in the physical (Hermitian) action defined in mo- 
mentum space. The incompatibility of these transforma- 
tions implies a change in the integration element of the 
partition function, which is given by 

II V Xi{k, l)V X 2(k, l)V X l(k, l)V x \{k, I) 
= niex P (W2)] 2 

LO.K 

x Y[V X i{K,uj)Vx2{K,cj)Vxi{K,cj)Vx\(K,oj) . (29) 



The integration element thus acquires an overall phase 
factor, which will drop out in the calculation of expec- 
tation values. This proves that expectation values cal- 
culated with the real-space action will be well defined 
and result in physical values. We attribute this to the 
fact that all steps in our derivation have a clear physical 
interpretation and motivation, and are not mere mathe- 
matical tricks. 



D. Adding gauge fields on the spacetime lattice 

We now continue this line of reasoning to formulate an 
action for chiral fermions interacting with gauge fields. 
Just as in the case of the Hamiltonian lattice, the way 
to go is by making the lattice for the gauge fields twice 
as fine as the lattice for the fermions. There will be four 
gauge field unit cells in one fermion unit cell. We label 
the parallel transporters living on the links as 



U x (k,l)=V 



U t (k,l)=V 



exp i / MxMt )^ 

;;+i/2)a t 
exp 1 1 J A t (ka x ,t)dt 

I a t 



, (30a) 



, (30b) 



with V being the path ordering operator, and this for k 
and I being either integers or half integers. 
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The correct way to add these transporters to the fermion action (19 1 is completely specified by definition (24 1. The 
result is given by 



S = a t Y,x\(k,l-l) 

k,l 



+xl(k,i-i) 



u t -\k, i - y a )xi(fc, i) - u t (k, i - i) X i(k, i - 1) 



a t 



. U-\k, I - i/2)x 2 (fc, I ~ 1) - U x (k - 1/2, 1 - y 2 )x 2 (fc - 1, / - 1) 



, u t -\k + y 2 , Q X2 (fc, Q - + y 2 , ; - i/2)x2(fc, j - i) 

a* 



(31) 



. u-\k + U)xi(k + i,0- Oxi(fe, 



The action in (31) still has the discrete symmetries A. Staggered fermions on the Hamiltonian lattice 



P, C, T and chiral symmetry, when accompanied by the 
appropriate transformation laws for the parallel trans- 
porters. In these transformation laws, it is important to 
include the possible shifts over half a distance a x and/or 
a t , as is required in the definitions of chiral transforma- 
tions and charge conjugations. These are again effects of 
order 0(a), which do not exist in the continuum limit. In 
the continuum theory for example, a chiral transforma- 
tion is only defined on the fermion fields and has no effect 
on the gauge fields. However, this discrepancy does not 
change the content of these transformations. Not only is 
the required shift of order 0(a) , but since it is a pure shift 
it has no influence on the the continuum theory which is 
required to have translational invariance. The only im- 
plication of these transformation laws is that chiral trans- 
formations and charge conjugations on the lattice can not 
be defined without performing such a shift. 



III. MINIMALLY DOUBLED, CHIRAL 
FERMIONS ON THE STAGGERED LATTICE IN 
(3 + 1) DIMENSIONS 



We are ready to apply the analogous approach to the 
chiral fermion action in (3+1) dimensions. We start again 
with the Kogut-Susskind construction on the Hamilto- 
nian lattice. However, fermion doublers will now already 
exist without discretizing time. We will briefly elabo- 
rate on the reason why we can not circumvent those by 
staggering the spinor components. Nextly, we will pro- 
ceed to the action formalism and define a unit cell in 
discrete space time which does not introduce extra time 
doublers. As such, we will have defined a non-Hermitian 
action which describes two chiral fermions. 

Since gauge fields can be added in a completely par- 
allel way, we will not repeat this argument and focus 
completely on the free, massless fermion field. 



We consider the three dimensional discretized space 
lattice. Keeping our approach to staggered lattices in 
mind, if the evolution equation for one field depends on 
the spatial derivative of another field, the second field's 
discretization position should be shifted over half a lattice 
distance in that spatial direction. The Hermiticity of the 
Dirac operator ensures that the evolution equation of the 
second field will depend on the same spatial derivative of 
the first field, which is required for the consistency of this 
approach. We choose to work with the standard represen- 
tation of the 7-matrices, although nearly any represen- 
tation is possible. Exceptions are those representations 
that result in the evolution equation of a field compo- 
nent being dependent on the spatial derivative of that 
same component. This is the case if one of the matrices 
a 1 = 7°7 l for i = 1,2,3 has non-zero diagonal entries, 
which is the case in i.e. the chiral representation or the 
Majorana reprensentation. 

In the standard representation, we start by placing tpi 
in the center of the unit cell. Since dipi jdt depends on 
dips/dz, ip3 should be discretized half a lattice distance 
a z further on in the z direction. The evolution equation 
of ipi also depends on the z-derivative as well as the y- 
derivative of ip4 . Continuing this line of reasoning results 
in the unit cell of FIG. [4] where all spinor components 
appear two times in the unit cell. We have the following 
rules for the discretization points: 

ipi : (ia>x,ja y , ka z ) 

((i + 1 /2)a x , (j + 1 /2)a y , ka z ) 
tp 2 : ((i + y2)a x ,ja v , (k + l /i)a z ) 

{ia x , (j + y^a-y, (k + 1 /2)a z ) 
ip 3 : {ia x ,ja y , (k + 1/2)02) 

((i + 1/2K, (j + y 2 )oy, (k + 1/2)0,) 
V> 4 : ((i+ 1 /2)a x ,ja v ,ka z ) 

fax, U + 1 h)a y , ka z ) 

Exactly the same construction was found by Susskind 
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in [5], again by introducing an auxiliary lattice which 
is twice as fine. Since we now have 8 independent fields 
within one unit cell, the dispersion relation in momentum 
space is found to be 



ui 2 — m 2 sin z 



(k x /2)- — sm 2 ( Ky /2) 



■ sin 2 (K z /2) 



(32) 



= 0. 



While there is only one low-energy region in the momen- 
tum space (k x ,k v ,k z ) £ {—n,+n{® 3 , it is fourfold de- 
generate. A factor two of this degeneracy is due to the 
different spin orientations. The other factor two points 
at two different fermion flavours, which is a logical con- 
sequence of having twice the number of fields in the unit 
cell. Susskind worked out the symmetries of this two 
fermion system in detail [20 and illustrated how the two 
independent fermion flavours could be reobtained. We 
will not repeat his argument. Let us only point out that 
a chiral transformation is obtained by shifting the lattice 
over half a lattice spacing a z in the z direction. This does 
indeed interchange the fields -01 V3 an d V"2 V'4- 



A; ct ; 



V>2 




V>3« 






m 

m 


m 













FIG. 4: Unit cell for the discretization of the fermion field 
components on the three-dimensional Hamiltonian lattice. 

Another way to look at the fermion doubling in this 
unit cell is by constructing the Voronoi cell or Wigner- 
Seitz cell, which is in fact the physical unit cell in which 
every field appears only once, and which is given in 
FIG. [5] For this smaller cell, the physical domain in 
momentum space is enlarged and a second low-energy 
region is discovered around the corners of the Brillioun 
zone, corresponding to the second fermion flavor. Both 
pictures are of course completely equivalent. 

Having this structure on the Hamiltonian lattice, we 
can proceed to the action on a discrete spacetime grid 
by performing the completely analoguous approach as in 
(1+1) dimensions. But let us first elaborate on the fact 
that the fermion flavor pair can not be reduced on this 
Hamiltonian lattice. The Dirac spinor in (3+1) dimen- 
sions has 4 components, just as the 4-tuple (p, v x , v y , v z ) 
in acoustics. However, the placement of the acoustics 
variables as given in FIG. [T] is not simply a trick to make 
the math work. This construction follows naturally from 
the physical interpretation of a scalar density p and a 



j a. 





+2tt 



(a) 



FIG. 5: Unit cell (dashed lines) and Voronoi or Wigner-Seitz 
cell (full lines) for the discretization of the fermion field com- 
ponents on the three dimensional Hamiltonian lattice, in up- 
per sight, (a) and the corresponding Brillouin zones in mo- 
mentum space (b). 



vector field v, and gives as such a clear representation 
of the conservation equation dp/dt + V • v = in a 
discrete setting. Such an interpretation is not possible 
for the four spinor components ipi, since spinors have no 
obvious classical, geometrical representation. 

We can try to find a representation of the 7-matrices 
for which a unit cell similar to that of the quadruple 
(p, v xi v y ,v z ) is useful. This would be the case if the time 
derivative of ip\ depends on the x-derivative of tp2, the 
y-derivative of V3 an d the z-derivative of V'4- Altough 
it is possible to construct three a-matrices with the re- 
quired properties, such as a 1 — \®a x , a 2 — o x ®a z and 
o? = <j x ® o y , it is now no longer possible to construct a 
matrix (3 which anticommutes with all a's. Also, trying 
to calculate 7 s ~ a 1 a 2 a 3 yields an answer which is pro- 
portional to the unit matrix. These conclusions clearly 
indicate that the fermion doubling problem should be 
completely attributed to the spinorial character of the 
field components and not, as is often thought, to the fact 
that the Dirac equation consists of a set of first order par- 
tial differential equations. First order partial differential 
equations introduce no problems whatsoever for acoustic 
or electromagnetic problems on space(time) lattices. 



B. Local lattice action for a minimally doubled 
pair of chiral fermions 

To construct a local action on a spacetime lattice, we 
should first look at a way to discretize time in the classi- 
cal equations of motion. Since the time derivatives of Vi 
and V2 depend only on (space derivatives of) V3 an d ip4, 
and vice versa, we should shift the discretization points of 
V3 and V4 m the time direction. In addition, the (1 + 1)- 
dimensional case has taught that the discretization points 
of the complex conjugate fields Vj should be shifted over 
half a lattice spacing at in the time direction. The re- 
sulting unit cell is depicted in FIG. [6] 
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The unit cell in FIG. [6] corresponds to the definitions 

Xla(hj,k,l) 



Xib(i,j,k,l 



yj Cl x (Xy(X z 
X2a(hj,k,l 

yj (l X (ly(X Z 

X2b(i,j,k,l 



yj Q, X CLy(l Z 

Xsa(i,j,k,l 



X3b(i,3,k,l 



yj Cl x (Xy(X z 

X4a(i,j,k,l 

yj (l X (ly(X Z 

X4b(i,j,k,l 



yJd~ X ~dy~Q~ Z ~ 



^i(t + ya,i + i/2,A5,0 

^(* + 1 /2,i,A + 1 /2»0 
i) 2 (i,j + 1/2, k+ 1/2,0 

V'3(i + 1 /2,J + 1 /2,fc+ 1 /2^ + 1 /2) 
V> 4 (« + 1 /2,J,fc^+ 1 /2) 



yjd~^dyd~z 
x\ b (i,j,k,l 



v/ Q'xQ>yQ>z 

Xlaihj'k, 1 
yj CL x (lyCL z 

x\ b {i,j,k,l 



X\a(hh k ,l 



yj Q x CLy(l z 

xl b (i,j,k,l 



yj Oj X CiyCL z 

x\ a {i,j,k,l 

yj Q j x® j y@ i z 

x\ b (i,j,k,l 



yJd~x~h~y~QT z 



3, k,l+ 1/2) 

^i(i + V2,i + V 2 »^^ + V 2 ) 

V'2(j+ 1 /2,J,fc+ 1 /2,/+ 1 /2) 
V'2(j,J + 1 /2,fc+ 1 /2^ + 1 /2) 
^(i,i,fc+l/2,/ + l) 

+ 1/2, J + Va,fc + i/2,i + l) 

^|(i+l/2,i,fc,Z + l) 

^(i,j + i/ 2 ,M+l) (33) 



in which the factor a^, a y , a z and a t in the arguments of the fields ipi were omitted for brevity. With these fields, we 
write down the discrete action 



S = \a t ^2x\a(h3'k,l) 

ijkl 



+xL(m>M) 



Xia(i,j,k,l+ 1) -Xia(i,3,k,l) , X3a(i,3, k, I) - X3a(i,3, k - 1, 



a t 



+xL(hJ\k,l-l) 



+Xla(h3\k,l- 1) 



Xia(i,3,k,l) - Xia(i - 1,3, k,l) .Xib(i,3,k,l) - Xib{i,j ~ 1,M) 



X2a(^ J, fc, Z + 1) - X2a(», J, fc, Z) _ Xia{hj,k + 1,1) - X4a(i , j , k , I) 

a t a z 
X3a(i + l,j, k,l) ~ X3a(i,3,k,l) | . XgbfejvM) - X3fc(»,j ^ lvM) 

X3a(i,3,k,l) - X3a(hj,k,l - 1) Xia(i,3,k + l,l) - Xia(i,3,k,l) 

at a z 
X2a(i, j, k, I) - X2a(i — l,j, k, Q _ { X2b(i, j, k, I) - X2b(i,j - l,fc,Z) 

CL X Qjy 

X4a(i,3,k,l) ~ X4a(i,3,k,l ~ 1) _ X2a(i,j,k,l) - X2a(i,j,k- 1,1) 

a t a z 
Xia(i + l,3,k,l) - Xia(i,3,k,l) | . Xu(i, j, k, I) - Xu(i, j - 1, k, I) 



(34) 



+ analogous expressions for the ^-fields. 



I 

Symmetry transformation which leave this action in- (l + l)-dimensional case and will not be written down ex- 
variant can easily be deduced from the continuum pre- plicitly. It is furthermore obvious that no new doublers 
scription and the definitions (33 1, along with a possible were introduced by this approach to time discretization, 
shift over half a lattice distance if the mapping of the We thus propose the action in (34 1 as a non-Hermitian 
fields would mismatch. A chiral transformation is ac- action for chiral fcrmions with ultralocal approximations 
complished by shifting the lattice over half a distance a z for the derivatives, which has the discrete symmetries P, 
in the z direction and half a distance a t in the t direc- C and T, and which describes only two fermion flavours, 
tion. All of this is a straightforward generalization of the 
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ia x ,ja y X >V 



FIG. 6: Unit cell for the discretization of the fermion field 
components on the spacetime lattice in (3+1) dimensions. 

IV. CONCLUSION 

We share the opinion of the authors of [17] that a 
Hermitian lattice action for minimally doubled, chiral 
fermions can not have all the required discrete symme- 
tries P, C and T. However, we argue that we can solve 
this problem by losing the Hermiticity condition at order 
0(a) for the discrete action, without this having to yield 
a physically ill-defined theory. We have illustrated a pos- 
sible approach by reinterpretating the staggered lattice 
approach in a language borrowed from the FDTD com- 
munity. In addition, we believe that these results imply 
a more stringent version of the Nielsen-Ninomiya theo- 



rem. A Hermitian lattice action for local, chiral fermions 
which has P, C and T symmetry will result in at least 
four fermion flavours. If the action is non-Hcrmitian, 
the Nielsen-Ninomiya theorem is no longer applicable, 
but, from our point of view, not all doublers can be re- 
moved without introducing non-physical effects or break- 
ing symmetries. The reason for this is that physicality 
requires the theory to be Hermitian for a t — > 0, and two 
different fermions will survive on the Hamiltonian lattice. 
This pair of fermions cannot be further reduced, due to 
the spinorial character of the field components. 

What we have given are not proofs of these statements, 
but some easy to understand arguments obtained by rein- 
terpreting this problem using the staggered lattice ma- 
chinery as it is applied in FDTD simulations. It remains 
to be seen wether these arguments are generally valid, 
and if these statements can be proven on a more rigorous 
level. 
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